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ABSTRACT 

A novel method to characterize the topology of the early-universe intergalactic medium during the 
epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is ana- 
lyzed through calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen-ionizing 
photons and ionization-rate equations are calculated in a suite of numerical simulations under various 
input parameters. The 2D genus is calculated from the mock 21-cm images of high-redshift universe. 
We construct the 2D genus curve by varying the threshold differential brightness temperature, and 
compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus 
curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate 
between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. 
Choosing the right beam shape of a radio antenna is found crucial for this analysis. Square Kilometer 
Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though 
some deterioration of the data for this purpose is unavoidable under the planned size of the antenna 
core. 

Subject headings: cosmology: theory - cosmology: observations - intergalactic medium - large-scale 
structure of universe - methods: statistical - radiative transfer 



1. INTRODUCTION 

The epoch of cosmic reionization (EOR) commences 
with the birth of first astrophysical, nonlinear objects 
such as the first stars and miniquasars. These sources 
of radiation create individual H II bubbles by emitting 
hydrogen- ionizing photons. This is then followed by the 
end of cosmic reionization when all percolating individual 
H II bubbles fully merge with one another and almost all 
hydrogen atoms in the universe become ionized. 

Direct observations of EOR have not yet been made, 
even though several indirect observations imply that 
(1) the end of reionizat ion occurred at a redshift z ~ 
6.5 (e.g. iFan et al.l [20Q6h . (2) the intergalactic medium 
remained at high te mperature during and after EOR 
(|Hui fc Haima n 2003) , and (3) this ep och started no later 
than z ~ 11 (jKomatsu et al I [2QQ9h . Theories suggest 
that the growth of H II bubbles proceeded inhomoge- 
neously (in patchy way) due to the clustered distribu- 
tion of radiation sources, and the global io nized frac- 
tion ( x) increased monotoni c ally in time (e.g. Illiev et al.l 
[20071 . iMcQuinn et "all I2007L IShin. Trac fcCei3 120081 : a 
few models predict, however, non-monotonic increase in 
the global ionized fraction (x) due to the possible re- 
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combination after regulated star formation, followed by 
eme rgence of stars of higher toleranc e to photoheating, 
e.g. ICen1l200l iWvithe fc Lo^l2003h . 

Observation of the redshifted 21-cm line from neutral 
hydrogen is one of the most promising methods for the 
direct detection of EOR, and the Cosmic Dark Ages that 
precedes this epoch as well. Both temporal and spatial 
fluctuations in the 21-cm signal are believed to be strong 
in general, thus easy to detect, during EOR due to the in- 
homogeneous growth of H II bubbles and relatively weak 
foregrounds at high frequencies. The strongest signal will 
come from the late dark ages, or very early EOR, when 
the spin temperature will be much lower than the cosmic 
microwave background radiation (CMB) due to the Ly- 
a coupling to yet unheated intergalactic medium (IGM), 
even though stronger foregrounds at lower frequencies 
will be an obstacle to observation. There are several 
large radio interferometer arrays which aim to detect 21- 
cm signal from EOR. These projects include Giant Me- 
terwave Radio Telescope^ (GMRT), Murchison Widefield 
Arra;y0 (MWA), Precision Array for Probing the Epoch 
of Reionizatiorf) ( PAPER: [Parsons et al.ll2010D . 21 Cen- 
timeter Arra£3 (21CMA; formerly known as PaST), the 
LOw Frequency ARra\f rT l (LOFAR), and Square Kilome- 
tre Arra£3 (SKA). 

The lack of direct observations of EOR results in poor 
constraints on the history of cosmic reionization. Theo- 
retical predictions for the history of cosmic reionization 
are made by semi- analytical calculations or fully numer- 
ical simulations. These studies select a variety of input 
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parameters, among which the most important one is the 
properties of sources of radiation. The mock 21-cm data 
produced by such studies will be compared to future ob- 
servations to constrain, for example, emissivity of high- 
redshift sources of radiation. 

The patchy, 3-dimensional 21-cm radiation background 
can be analyzed in various ways such as the 3D power 
spectrum, 2D power spectrum, distribution of H II bub- 
ble size, cross- correlation of ion i zed fraction and o verden- 
sity, etc (e.g. Illiev et al.l [20071 : IZahn et al-l l2006h . Each 
analysis method adds to the capability to discriminate 
between different reionization scenarios. They are usu- 
ally complementary to each other, allowing one to un- 
derstand the underlying physics better the more tools 
are available. Therefore, it is always favorable to have as 
many different tools for data analysis as possible. 

Here we introduce a novel approach to the analysis of 
21-cm radiation background from EOR. We characterize 
the geometrical property of the distribution of neutral 
IGM and H II bubbles by using the topology of 21-cm 
differential brightness temperature field. For this pur- 
pose we measure the 2D genus of a series of snapshots 
of the high-redshift universe that are predicted in dif- 
ferent models of reionization. By varying the threshold 
value for the differential brightness temperature, we con- 
struct 2D genus curves at different redshifts under differ- 
ent reionization models. These models are simulated by 
self-consistent calculation of radiative transfer and rate 
equations over our simulation box. 

Recently, similar methods for studying topology of 
high-r edshift IGM have been su ggested by IGleser et all 
(|2Q06f ). iLee et al.l (|2008f ). and iFriedrich etHI (|2Q10f ). 
They calculated the 3D genus using eith er the neu- 
tral (IGleser et al.l |2QQfl ILee et al.l l2QQ8h or ionized 
(|Friedrich et al.l l2010f ) fraction of IGM. These studies 
show that the 3D genus of the underlying neutral (or ion- 
ized) fraction reflects the evolutionary stages of cosmic 
reionization. The 3D genus is also found to be useful 
in discriminating between different reionization scenar- 
ios. Albeit the similarity of our work to these studies, 
there are several fundamental differences. First, we cal- 
culate the 2D genus instead of the 3D genus. Second, we 
use an observable quantity, the 21-cm differential bright- 
ness temperature, which is more directly applicable to 
real data observed by future radio antennae, when cal- 
culating the 2D genus. Third, we explicitly seek for the 
possibility to apply our method to the data to be ob- 
served by SKA, with appropriate choice on the filtering 
scales for the 21-cm differential brightness temperature. 

We emplo y the 2D genus method developed by 
iMelott et all (|1989D . The method of genus topol- 
ogy was originally designed to test mainly for 
the Gaussian random-phase nature of the primor- 
dial density field in 3D (fGott. Melott fc Dickinson! 

19861: Hamilton, Gott fe Weinberg! 119861: 

Gott . Weinberg fe Melottl I1987D or in 2D (jMelott et all 
1989). The 2D case has been studied fo r a variety of cos- 
mological data sets: on redshift slices dPark et al.l 119921 : 
lCollevlll997l : iHovle. Vogelev fe Got t 2002), on sky maps 
([Gott et al. 1992; Par k. Gott fe Choi 2001) on the CMB 
iCole s 1989; Got t et aLNl 990: Sm oot et al . 1992; Kogutl 
[1991; Ep gut et al. 1996; iCollev. Gott. fe Parkl [1991 
iPark et alT [l998. 200 1]), and on the HI gas distribution 
in galaxy disks ([Kim fe Par Id 120071 : iChepurnov et al.l 



EooD. 

Because 2D genus is strongly affected by the number 
and distribution of H II bubbles as well as the distribu- 
tion of underlying density field, our method may provide 
an insight about how reionization proceeded and was af- 
fected by properties of sources. We also seek for the 
possibility to apply our method to observed data from 
SKA, which so far has the highest sensitivity proposed. 

The plan of our paper is as follows. In Section [2] we 
describe our numerical simulations and the basic proce- 
dure to calculate the 21-cm radiation background. In 
Section [3] we describe how the 2D genus is calculated. In 
Section [U our simulations results are analyzed and possi- 
bility of using SKA observation results are investigated. 
We summarize our result and discuss future prospects in 
Section \5\ 

2. SIMULATIONS OF REDSHIFTED 21-cm SIGNALS 
2.1. N-body simulations 

We completed a 2048 3 particle simulation in a con- 
cord ance ACDM model with WMAP 5-year parame- 
ters (jDunklev et al.l 12009): Q m = 0.258, ^ A = 0.742, 
Q b 0.044, n s 0.96, a 8 0.79, and h 0.719, where 
^m,^A,^6 are the density parameters due to matter, 
cosmological constant, and baryon, respectively. Here, 
n s is the slope of the Harrison- Zeldovich power spec- 
trum, and o~ 8 is the root-mean-square (rms) fluctuation 
of the density field smoothed at 8/i _1 Mpc scale. The 
cubic box size is 66/i _1 Mpc in a side length. The sim- 
ulation uses GOTPM a hybrid PM+T ree AT-body code 
(jDubinski et al.l 120041 : iKim et al.l 12009). The initial per- 
turbation is generated with a random Gaussian distribu- 
tion on a 2048 3 mesh at Z{ = 500. The force resolution 
scale is set to f e = 3.2/i _1 kpc in the comoving scale. 
The total of 12,000 snapshots, uniformly spaced in the 
scale factor, are created. The time step is predetermined 
so that the maximum particle displacement in each time 
step is less than the force resolution scale, f e . 

We then extract friend-of-friend (FoF) halos at 84 
epochs: 24 epochs with five million year interval from 
z = 40 to 20, and 60 epochs with ten million year sep- 
aration between z = 20 and z = 7. The connection 
length is set to 0.2 times the mean particle separation to 
identify cosmological halos. Using this method, we find 
all halos with mass above 1O 8 M (corresponding to 30 
particles). At each epoch, we calculate matter density 
fields on 2048 3 mesh using the Triangu lar- Shaped- Cloud 
(TSC) scheme (Efst atFiou et al.l[l9 85). As the radiative 
transfer through the IGM does not need to be run at too 
high resolution, we further bin down the original density 
fields to 256 3 after subtracting the contribution of col- 
lapsed dark halos. We transform this dark matter den- 
sity to baryonic density by assuming that in the IGM 
baryons follow the dark matter with the mean cosmic 
abundance. From the list of collapsed halos, we form a 
"source catalog" by recording the total mass of low-mass 
(M < 1O 9 M ) and high-mass (M > 1O 9 M ) halos in 
those radiative transfer cells containing halos. 

We also ran a similar A^-body simulation with the same 
cosmological parameters but a different configuration of 
512 3 particles on a 512 3 mesh on a 64 h~ x Mpc box using 
GOTPM, resulting in somewhat poorer halo-mass reso- 
lution of 2 x 10 9 Mq. We created a binned-down density 
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Table 1 

Simulation parameters. Cases 1 and 2 have only high-mass sources. Cases 3 and 4 have both high-mass and low-mass sources. Numbers in 
the "Model" are / 7 's for high-mass halos, and if any follows, are / 7 's for low-mass halos (with "S" representing "self-regulated"). 
Corresponding g 7 's are also listed. Radiative transfer is calculated on the number of meshes listed under "RT mesh." Note that even 
though case 1 has the same g 7 as cases 3 and 4, its minimum mass M m i n does not reach the low-end (10 9 Mq) of the other two, thus 

making its reionization history end much later. 



Case Model L box (Mpc/h) TV-body mesh RT mesh M min (M©) / 7 , h igh / 7> iow At (Myr) z 50% z 1% 

1 f250 (gl25) 64 512 3 128 3 2 x 10 9 250 - 20 9.7 11.5 

2 f5 (g2.5) ' 64 512 3 128 3 2 x 10 9 5 - 20 7.1 9.7 

3 fl25_1000S (gl25_1000S) 66 2048 3 256 3 10 8 125 1000 10 12.0 17.3 

4 fl25_125S (gl25_125S) 66 2048 3 256 3 10 8 125 125 10 11.6 15.2 



field on a 128 3 mesh, similarly, on which radiative trans- 
fer is calculated. 

We use the results from the high resolution (2048 3 TV- 
body) run for two cosmic reionization simulations (Case 
3 and 4), and those from the low resolution (512 3 TV- 
body) run for the other two (Case 1 and 2). For each 
cosmic reionization simulation, we choose a fixed time 
step At. These parameters are listed in Table [TJ 

2.2. Reionization simulations 

We performed a suite of cosmic reionization simula- 
tions based upon the density field of IGM and the source 
catalog compiled from the N-body simulations, described 
in Section EEJ We used C 2 -Ray (jMellema et all I2005D 
to calculate the radiative transfer of H-ionizing photons 
from each source of radiation to all points in the simula- 
tion box and the change of ionization rate at each point, 
simultaneously. The spatial resolution of the binned- 
down density field is the actual radiative-transfer reso- 
lution, which is described in Table [T] for each case. At 
each TV-body simulation step, we produce corresponding 
3-dimensional map of ionized fractions. 

We vary physical properties of sources by changing the 
parameter f 1 = / esc f* N{ over simulation, where f esc is 
the escape fraction of ionizing photons, f* is the star 
formation efficiency, and Ni is the number of ionizing 
phot ons emitted per s tellar baryon during its lifetime At 
(e.g. Illiev et a l. 2007). / 7 determines the source prop- 
erty in such a way that / 7 M tot photons are emitted from 
each grid cell during At, where M to t is the mass of halos 
inside the cell. At is also used as the time-step for finite- 
differencing radiative transfer and rate equations. If At 
varies among simulations, fixed f 1 will generate different 
numbers of photons accumulated. Therefore, it is some- 
times preferable to use a quantity som ewhat blind to 
At. q 7 = / 7 /(At/10Myr), as defined in iFriedrich et al.l 



(|2010f h is such a quantity, which is basically the emissiv- 
ity of halos. 

We also adopt two different types of halos: low-mass 
(10 8 < M/M© < 10 9 ) and high-mass (M/M© > 10 9 ). 
Both types are "atomically-cooling" halos, where colli- 
sional excitation of atomic hydrogen Lya line predomi- 
nantly initiates cooling of baryons, reaching the atomic- 
cooling temperature limit of ~ 8000 K, and is followed 
subsequently by H2-cooling into protostellar cloud. One 
main difference exists in their vulnerability to external 
radiation: low-mass halos, when exposed to ionizing ra- 
diation, the Jeans mass after photoheating overshoots 
their virial mass such that the gas collapse is prohib- 
ited, while high-mass halos have mass large enough to 



be unaffected by such photoheating. We thus adopt a 
simple "self-regulation" criterion: when a grid cell ob- 
tains x > 0.1, we fully suppres s star formation inside 
low-mass halos in the cell (e.g. Illiev eFaII[2007l ). We 
also use in general different / 7 's for different types, when 
we include low-mass halos. Cases 1 and 2 did not allow 
us to implement low-mass halos due to the limited mass 
resolution, while cases 3 and 4 did, thanks to the high 
mass resolution. 

Finally, we note that halos in even smaller mass range, 
10 5 < M/M© < 10 8 , or minihalos, are not consid- 
ered in this paper. This can be justified because their 
impact on cosmic reionization, e specially at late stage, 
are believed to be negligible ( e.g. H aiman, Abel fe R ees 
120001 : lHaiman fc Brvanl I2006D . Nonetheless, the earli- 
est stage of cosmic reionization must have started with 
minihalo sources, most probably regulated by inhomo- 
geneous H2 Lym an- Wer ner (LW) radiat ion background 
(jAhn et al.lfcooa see also lDiikstra et al.ll2008l for the im- 
pact of LW background on growth of supermassive black 
holes). A preliminary but self-consistent simulation of 
cosmic reionization including minihalos and LW feed- 
back has been achieved very recently, showing that about 
~ 10% of global ionization of universe can be completed 
by minihalo sources at high redshifts (|Ahn et al.l 120101 : 
Ahn et al. in preparation). 

2.3. Calculation of 21-cm differential brightness 
temperature 

The hyperfine splitting of the ground state of hydrogen 
leads to a transition with excitation temperature X* = 
0.068 K. The relative population of upper (ni) and lower 
(no) states are determined by the spin temperature T5, 
such that 

Several pumping mechanisms determine Ts as follows: 



^CMB 



VcT, 



l + y a + y c 



(2) 



where Xcmb, T ai Tk, y a , and y c are the cosmic mi- 
crowave background (CMB) temperature, Lya color tem- 
perature, the kinetic temperature of gas, the Lya cou- 
pling const ant, and the collisional coupling c onstant, re- 
spectively (jPurcell fc Fieldl [1956L [Fiildl fl959l) . 

21-cm radiation is observed in emission against CMB 
if Ts > Xcmb, or in absorption if Ts < Xcmb- This is 
quantified by the 21-cm differential brightness tempera- 
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ture (e.g. IMorales fc Hewittll200l 

6T b = Ts - rcMB (1 - e -) , (3) 

J- \ z 

if 21-cm lines from neutral hydrogen at redshift z is trans- 
ferr ed toward us. T he optical depth r is given by (see 
e.g. Illiev et al.ll2002l for details) 

x fAWikV"\ 1+il . (4) 

becomes appreciable enough to be detected easily 
only when Ts > Tcmb ot Ts <C Tomb- The former 
limit is believed to be reached quite early in the his- 
tory of cosmic reionization, because X-ray sources can 
easily heat up the gas over cosmological scales with 
small optical depth, su ch that Ts ~ Tk ^> Tcmb (e.g. 
iCiardi fe Madaul l2QQ3h . There exists another limiting 
case where Ts is determined only by collisional pumping 
(y a C 1), which is usually the cas e in the Cosmi c Dark 
Age s as studied by, f or instance, iShapiro et al.l (j2QQ6h 
and|K 1m fc Ped(|2010h . In this paper, we simply assume 
early heating of IGM such that Ts ^> Tcmb- Equa- 
tion (j3]) and (|4j) then give the limited form of JT5, 

We simply generate 3D maps of ST^ on the radiative 
transfer grid over all redshifts using Equation ([5j). Note 
that in the limit Ts ^> Tcmb, ^T^ oc (1 + S)(l — x) 
at a given redshift. This limit allows one to detect ST^ 
relatively more easily than the case where Ts ~ Tcmb, 
and generates a 3D map proportional to the underlying 
density field as long as gas is fully neutral (x = 0). Of 
course, patchy evolution of H II bubbles will be strongly 
imprinted on this ST^ map as well. 

3. 2D GENUS OF 21-cm DIFFERENTIAL BRIGHTNESS 
TEMPERATURE 

In 3D, the genus of a single connected surface is iden- 
tical to the number of handles it has. In 2D, if a field of 
a variable T is given, the 2D genus for a threshold value, 
T th , is given by 

£2 D (T th ) = JV+(T th ) - M(Tth) , (6) 

where 7V + (T t h) and 7V_(T t h) are the number of connected 
regions with T > T t h (hot spots) and T < T t h (cold 
spots), respectively. In general, as one varies T t h, #2D 
changes. Therefore, #2D is a function of Tth, and its 
functional form is called the 2D genus curve. In our 
case, we calculate the 2D genus curve of the field of 21- 
cm differential brightness temperature JT5, by varying 
the threshold value 5Tb, th- 

The 2D genus curve of a field can work as a character- 
istic of the field. A useful template for this is a Gaussian 
random field, because its 2D genus is given analytically 



by the relation 

#2D, Gauss OC ^exp(-Z/ 2 / 2 ) , (7) 

where 

„= T *-W (8) 

is the deviation from the average (T) in units of standard 
deviation o~t- One can then compare the 2D genus curve 
of a given field to #2D, Gauss to see how close the field is 
to a Gaussian random field, for example. 

The 2D genus is also equal to the integral of the curva- 
ture along the contours of T = T t h, divided by 2tt. This 
is because when a curvature is integrated along a closed 
contour around a hot spot (cold spot), its value is 2tt 
(— 2ir). When a 2D field is pixelized, a contour of T = T t h 
is composed of a series of line segments with turns occur- 
ring at vertices. On our rectangular 2D grid, every vertex 
shares four pixels (except for vertices on the edge and the 
corner). Under these conditions, the CONTOUR2D pro- 
gram is able to calculate the 2D genus by counting the 
turning of a cont our observed at eac h vertex of four pix- 
els in an image (jMelott et al.|[l989l ): 1/4 is contributed 
to the total genus from each vertex with 1 hot pixel and 
3 cold pixels, —1/4 from each vertex with 3 hot pixel and 
1 cold pixels, and zero otherwise. We do not calculate 
the modified 2D genus, #2D,eff = #2D — 2/, where / is the 
areal fraction of hot spots on the sky. ,q f 2D r eff is appro- 
priate when applied to a full-sk y field (|Gott et al.lll990l : 
iCollev fc Go"ttll2QQl iGott et aLll2QQ7h . while #2D is more 
appropriate for a relatively small, restricted region on the 
sky. Our simulation box is large enough to capture the 
typical size of H II bubbles (~ 10 — 20 Mpc comoving) 
during most phases of the cosmic reionization, but still 
too small to cover the whole sky and capture fluctuations 
in much larger scales. 

Our aim is to obtain 2D genus curve from the 
frequency- and angle- averaged differential brightness 
temperature at given frequencies. In order to compare 
different reionization models, we will choose those fre- 
quencies such that different reionization simulations yield 
the same global ionized fraction. We will then vary the 
threshold (5Tb, th to construct 2D genus curve ^2D(^Tb,th) 
for each case. 

The amplitude of 2D genus is roughly proportional to 
the field of view at a fixed redshift. Throughout this 
paper, the 2D genus we present is normalized to the field 
of view of size (1 degree) 2 . 

4. RESULTS 

4.1. Mock 21-cm sky maps 

Radio signals, observed by radio antennae, should be 
integrated over the angle, frequency and time in order 
to achieve the sensitivity required for specific scientific 
goals. Observing signals from high redshifts are chal- 
lenging, especially, because these signals are very weak 
while the required angle- and frequency-resolutions are 
relatively high. Therefore, very intensive integration in 
time (~ 1000 hours) is required for EOR observations in 
general. 

We therefore choose the beam size and the frequency 
bandwidth which we find adequate for generating dis- 
tinctive 2D genus curves from different reionization see- 
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Figure 1. 21-cm maps at redshifts z = 15.645, 12.184 and 11.325 for fl25_1000S. Top row: Slices which are frequency-averaged into one 
band of Av = 0.2 MHz, showing ionized regions (yellow) and neutral regions (blue), superimposed on the density field (light and dark for 
high-density and low-density regions, respectively). Middle row: Corresponding 2D 21-cm differential brightness temperature, averaged by 
a Gaussian beam of AO = V . Bottom row: The same as middle, but with a compensated Gaussian beam of the same AO. 

narios, under the assumption of very long time integra- bandwidth 
tion (~ 1000 hours or more). We generated mock 21- 
cm emission sky maps for our four reionization scenar- 
ios, f5, f250, fl25_125S and fl25_1000S, averaging ST b on 
our computational grid over the angle and the frequency. 
First, note that our box size corresponds to the frequency 



Az7box ^TTW 2 

-(4.10 MHz) ( ) ( 

V } ymh-iMpcj Vo.7i9y 



V2 



0.258 ) (, 15 ' (9) 
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and the transverse angle 



A#box = 



7 box 



7 box 



(l + z)D A (z) r(z) 

^box 



: (30.25 7 ) 



66 ft- 1 Mpc 



/10563.5 - 11863.5(1 + ^T 1/2 \ 
X V 10563.5 - 11863.5(15)- 1 /2 J ' ( > 

where v§ = 1420 MHz is the rest frame frequency of the 
line, Da(z) is the angular diameter distance, Lbox * s ^ ne 
comoving length of our simulation box, c is the speed of 
light and Hq = 100 /ikms -1 Mpc -1 is the Hubble param- 
eter at present. The numerator in the last parenthesis 
in equation (pTU]) is a fitting formula for the line-of-sight 
comoving distance r(z) in units of Mpc, under the cos- 
mological parameters used in this paper. 

We first calculate bare differential brightness tempera- 
ture on our computational grid by using Equation ([5j). 
We then integrated them over frequency with chosen 
Az/s. We also consider the Doppler shift due to peculiar 
velocity. Along the line-of-sight (LOS), the frequency- 
averaged 5Tb then becomes 



(5T b ) Au 



(11) 



where N Av = AL(Av) /l ce \\ is the number of cells corre- 
sponding to Av at a given redshift, fi is the fraction of 
a cell i entering the band after the Doppler shift, 



AL(Au) ■: 



Avc(l + zf 



' Z/ #0\AMl + ^) 3 + ^A 

: (3.22 Mpc) (-^=\ f^P _ 
v F J V 0.2 MHz J V 0.719, 



,0.258 ) 



-1/2 / 1 + 2 XV2 
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(12) 



is the comoving length along LOS corresponding to Az/, 
and / ce n is the comoving length of a unit cell. For a sta- 
tistical purpose, we generate four 21-cm emission maps 
for each LOS direction from a single simulation box by 
periodically shifting the box by a quarter of box size 
(^box/4). These maps will be almost mutually indepen- 
dent, because Lbox/4 ^ AL(Ais). We also take three 
different directions for LOS (a>, and z-direction) by 
rotating the box, such that we have in total 12 (almost) 
independent 21-cm sky maps from a single simulation 
box at a given redshift. 

Because the shape of the beam varies over different an- 
tennae, we choose two different types of beams to average 
the signal in angle: Gaussian and compensated Gaussian. 
We first convolve the 21-cm emission map from simula- 
tions with a Gaussian filter 



w Q (p) 



2tt<7 2 6XP V 2a 2 



(13) 



with a full width at half maximum (FWHM) AO = 
2oV2 ln2. We also use a compensated Gaussian filter 



given by 

Wcg(9) = 



27TO- 2 



^-^) eXP ( 1_ ^> 



(14) 



with a FWHM given by 



AO = 2(7^2(1 - LambertW(e/2)) (15) 

where LambertW(e/2) w 0.685. 

A Gaussian filter is widely used in literature to ap- 
proximate actual beams. In general, an interferometer 
will not retain the average (large scale) signal, but only 
record the fluctuations between minimum and maximum 
angular scales (depending on the largest and smallest 
baselines). As the largest scales to which the upcoming 
interferometers are sensitive are larger than our simu- 
lated images, a reasonable approach would be to sub- 
tract the mean signal from our images. However, for the 
determination of the genus this does not make a differ- 
ence other than shifting the threshold value, so we retain 
the average signal here. A compensated Gaussian filter 
roug hly mimics the beam of a compact interferometer 
(e.g. Mellem a~eTaLl 120061 ). even though real beams have 
sometimes much more complicated shape depending on 
the actual configuration of the antenna. We just take 
Wqg as an extreme case for "dirty" beams, clearly dis- 
tinct from Wq, because this filter is somewhat patholog- 
ical as it suppresses both small and large scale features. 
When > \J2cj, Wcg(#) < such that a data field con- 
volved with Wgg(&) ma y change its sign from point to 
point. 21-cm signals filtered this way would show seem- 
ingly unnatural H II bubble feature, or in a worse case, 
even absorption signals (ST^ < 0), even though Equa- 
tion ([5]) implies that 6T b > when T s > T C mb- 

Figure [T] shows how the actual 2D 5T b field (top rows) 
will be observed in different filtering schemes (middle: 
Wg; bottom: Wcg)- Three different epochs were chosen 
to represent early (volume weighted ionized fraction x v = 
0.04), middle (x v = 0.4), and late (x v = 0.9) stages. 
Wq generates reasonably filtered ST^ maps at all three 
epochs. At very early (x v = 0.04) and late (x w = 0.9) 
stages, fluctuation in 5T b is dominated by 1 + 5 and 1 — x, 
respectively, while at the middle [x v = 0.4), both 1 + 6 
and 1 — x contribute to the fluctuation. 

The compensated Gaussian beam adds ripple struc- 
tures to the emission signals. When a uniform field con- 
fined within a certain boundary is convolved with Wgg , 
this ripple becomes visible only along the boundary. As 
the ionized regions grow, therefore, a compensated Gaus- 
sian beam produces negative through of width similar to 
the beam FWHM at the boundaries of bubbles. In gen- 
eral, Wqg changes the topology of the 21-cm signal. For 
most of our analysis, therefore, we just use Wq. Fur- 
ther comparison of these two filters will be made in Sec- 
tion [4321 

4.2. Evolution, sensitivity, bubble distribution, and 
power spectrum 

The evolution of global volume-weighted ionized frac- 
tion x v varies significantly over different reionization sce- 
narios, as seen in Figure O First, the end of reionization 
in f5 and f250 cases occurs relatively later than fl25_125S 
and fl25_1000S cases, because these are tied only to high- 
mass (M > 2 x 1O 9 M ) halos, which collapse much later 
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Figure 2. Left: Evolution of the mean volume ionized fraction as a function of redshift. Middle: The mean differential brightness 
temperature with V FWHM Gaussian beam and 0.2 MHz bandwidth as a function of the mean volume ionized fraction. Right: rms 
fluctuation of the differential brightness temperatures filtered by the same beam and bandwidth as a function of the mean volume ionized 
fraction. 



than low-mass (10 8 < M/M© < 10 9 ) halos. Second, f5 
and f250 also show steeper evolution in x v , because there 
is no chance for slow evolution of total luminosity from 
self-regulation of small-mass halos. 

It is in general a better practice to compare different 
models at a fixed global ionized fraction than at a fixed 
redshift, partly because there are still too much freedom 
in the exact epoch of the end of reionization, the time- 
integrated optical depth to Thompson scattering of CMB 
photons, etc. More importantly, a fixed global ionized 
fraction among different models is reached by roughly the 
same number of H-ionizing photons accumulated from 
the beginning of EOR. We therefore use x v as the time 
indicator of the global evolution throughout this paper. 

Using x v as the time indicator, we first observe simi- 
lar trends in the evolution of 5Tb among different mod- 
els. The mean differential brightness temperature (5Tb) 
starts from (5Tb) « 30 — 40 mK and gradually decreases 
in time, as seen in Figure O The root-mean-square 
(rms) fluctuation of 5Tb, 5T& jrms , reaches the maximum 
(~ 8-11 mK) at x v « 0.4 - 0.6 (Figure E}, when the 
signal is filtered with Au = 0.2 MHz and AO = V . 

Albeit the similar trends, detailed analysis will be 
needed to discriminate between different models and ul- 
timately probe the properties of sources. The probability 
distribution function (pdf) of bubble size is a useful tool 
to understand the impact of these source properties. The 
bubble size is associated with the mass spectrum of halos 
and the ionizing efficiency / 7 , simply because it roughly 
reflects the total number of H-ionizing photons emitted 
into the bubble. In order to find the pdf of bubble size, 
we use a hybrid method. Fi rst, the size of a bu bble is de- 
termined by the method in IZahn et al.l ([20061 ): the bub- 
ble size is the maximum radius of a sphere from each 
simulation cell inside which the ionized fraction is over 
90%. This way, every cell in the box is associated with a 
bubble wit h a certain size. We then use the void-finding 
method by iHovle fc Vogelevl (|2002[ h which has been ex- 
tensively used for finding cosmological voids which are 
mutually exclusive in space. The bubbles are sorted in 
size from the largest to the smallest. The largest bubble 
is considered an isolated one. We then move to the next 
largest bubble and check if there is any overlap in vol- 
ume with the fist one. We iterate this over the bubble 
list until the overlap becomes less than 10%, which then 



registers another isolated bubble. This process is further 
iterated and we obtain the full list of mutually exclusive 
bubbles (to the extent of 10% overlap). 

In Figure 03 we show the pdf of bubble size obtained 
this way. We can observe clear distinction of f5 and f250 
from fl25_125S and fl25_1000S. In f5 and f250, a given 
ionized fraction is reached only by high- mass halos. For 
example, to reach x v ~ 0.6, the evolutionary stage of 
cosmological structure formation should enter a highly 
nonlinear phase because there is only one species (high- 
mass halos) responsible for such high x v . In contrast, 
in fl25_125S and fl25_1000S, both of the two different 
species (high- and low- mass halos) contribute to, for ex- 
ample, x v ~ 0.6. Therefore, the structure formation will 
be in less nonlinear stage than the former cases. Corre- 
spondingly, we can expect stronger clustering of sources 
for f5 and f250, and weaker for fl25_125S and fl25_1000S. 
Correspondingly, merger of bubbles would be stronger for 
the former and weaker for the latter, respectively. This 
is finally reflected in the relative contribution to a given 
x v from largest bubbles, as seen in Figure 03 Large bub- 
bles dominate x v in f5 and f250, while there are relatively 
smaller contribution to x v in fl25_125S and fl25_1000S. 

We also obtain 2D angular power spectra of the neutral 
fraction, 1 — x, for further analysis. This also provides 
an insight on the clustering scale of sources of radiation, 
or H II bubbles, through the location of peaks in the 
power spectrum. We follow the convention used in CMB 
analysis and constructed [£(£ + l)C^/27r] 1 / 2 in spherical 
harmonics, where £ roughly corresponds to 2tt/0. See 
Figure [H At x v = 0.04, it is clearly seen how cosmic 
reionization commences in each scenario. f5 and f250 
show peaks in £{£ + A)Cg at £ ~ 10 4 , or the comoving 
length scale x ~ 6 - 7 h~ l Mpc, while fl25_1000S shows 
a peak at £ ~ 3 x 10 4 , or x ~ 2 h' 1 Mpc, and fl25_125S 
shows a long uniform tail from 3 x 10 4 to ^ 4 x 
10 4 . The behavior of 2D angular power spectra indicates 
stronger merger of (otherwise) individual bubbles into 
larger scale at x ~ 6-7 h' 1 Mpc in f5, f250 and fl25_125S 
than in fl25_1000S. Even though fl25_125S has a tail 
up to £ - 4 x 10 4 (x ~ 1 ft -1 Mpc), the smaller / 7 (= 
125) in small halos than that of fl25_1000S (f y = 1000) 
makes these small halos less efficient. This is true at 
all redshifts: the power spectrum of fl25_125S is hardly 
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Figure 3. Left: The probability distribution functions of bubble size . We determined t he siz e of a bubble using the method in Zahn et al. 
(2006), and select mutually exclusive bubbles using the method in Hovle & Vogelev (2002). Right: The volume- weighted probability 
distribution functions of bubble size. 
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Figure 4. 2D angular power spectra of neutral fraction 1 — x. 

distinguishable from those of f5 and f250 any time except 
for higher power in the smallest length scale (highest £). 
In contrast, a distinctive clustering scale, £ ~ 3x 10 4 (x ~ 
2/*- 1 Mpc), is shown at all redshifts in fl25_1000S, due 
to the clustering of small halos of the highest efficiency. 

In short, both bubble pdf and 2D angular power spec- 
trum of neutral fraction may work as useful tools for un- 
derstandi ng th e nature of different reionization scenarios. 
In Section 331 when we analyze 2D genus curves of differ- 
ent cases, we will address these properties to see whether 
2D genus curves created from possible observations can 
bear similar fundamental nature as well. 

4.3. Genus properties 

We generate g^D^T^) from selected epochs at which 
x v = 0.04, 0.4, 0.6 and 0.9 as described in Section [3l 
The base 2D field of ST^ is averaged in frequency with 
Ai/ = 0.2, 1 and 2 MHz, and in angle with AO = 1', 2' 
and 3'. We vary Av and AO to quantify the competing 
effects of changing resolution and sensitivity, and also 
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Figure 5. 2D genus denned by the threshold differential bright- 
ness temperature within ~ 30' x 30 7 square (corresponding to our 
box size) in the sky assuming all the gas was neutral, filtered with 
V FWHM Gaussian beam and 0.2 MHz bandwidth (red). In com- 
parison, we also show the 2D genus of the differential brightness 
temperature with the same condition but from one-cell thick slice 
in our simulation (blue). S^D^^n) of the bare A/-body field is 
thus subject to smoothing over the size of the cell. This result is 
based upon the high-resolution N-body simulation (used for cases 
3 and 4). 



to make comparison with characteristics of future radio 
antennae. In angle-averaging, only Wq is used in all 
cases, except for fl25_1000S to which Wqg is also applied 
to understand the impact of beam shape on g 2 D- 

4.3.1. 2D genus of density fluctuation 

Is there any useful template genus curve to be com- 
pared with g2F>(5Ti))? There is such one indeed, which is 
#2D of an artificial field of ST^ which reflects the density 
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fluctuation only. Let us denote this quantity by 
ST btn = 5T b (x = 0) 

Distribution of ST^n will be Gaussian in the linear regime 
(\5\ <C 1), which will provide well-defined, analytical 
02D,Gauss(#Zb,n). Even when g2D(5T bjn ) starts to devi- 
ate from #2D, Gauss (<^6,n) due to nonlinear evolution of 
high density peaks, #2D(<^6,n) will be used to indicate 
overdense (5 > 0) and underdense (5 < 0) regions. Fig- 
ure [5] implies that ST^n > corresponds to 5 > 0, while 
5T b , n < to 5 < 0. 

g2T>(STi Jjn ), plotted in Figure El show how matter den- 
sity evolves. At z ~ 10 — 15, there already exist ex- 
tended tails into high 5Tb, ni because these correspond 
to high-density peaks that have evolved into nonlinear 
regime. Note that purely Gaussian random field gener- 
ates #2D, Gauss > which is symmetrical around v = 0, while 
high-density peaks (on high v) that evolved nonlinearly 
deviates from the Gaussian distribution to also make #2D 
deviate from #2D, Gauss- Because we will use filtered 5Tb 
maps, g2B of filtered 5Tb jU will be the template to be used 
throughout this paper. Filtering makes both the ampli- 
tude and width of g2T>(5Tb, n ) shrink (Figure [5]). Nev- 
ertheless, even after filtering, regions with positive #2D 
correspond to overdense regions and negative #2D to un- 
derdense regions, because the mean-density point (where 
the #2D(<^6,n) curve crosses the 5Tb axis) is almost un- 
shifted. This will work as a perfect indicator, if compar- 
ison is made with #2d(^1>), about how different the 5Tb 
field is from the underlying density field. #2D(#?6,n) can 
be obtained accurately for any redshift as long as cosmo- 
logical parameters are known to high accuracy, which is 
the case in this era of precision cosmology. 

4.3.2. Evolution of #2D and impact of beam shape 

We found that the evolution of g2D(5Tb) contained 
both generic and model-dependent features. We describe 
the generic features based on fl25_1000S case. We wil l 
describe the model-dependent features in Section 14.3.31 

We find that the change of g2T>(5Tb) clearly shows how 
reionization proceeds in time. The left panel of Figure [6] 
shows Wb-filtered g2D(5Tb, n ) of the underlying matter 
density field (blue dashed line) and g2u(5Tb) of the dif- 
ferential brightness temperature field of HI gas (red line). 
In the earliest phase of reionization (x v = 0.04), g2T>(5Tb) 
is much smaller than g2r>(5Tb, n ) at high temperature 
thresholds (5Tb, n > 50 mK). It is because the high- 
est density peaks have been ionized by the sources in- 
side them, and are dropped out of the neutral gas dis- 
tribution. It can be also noted that the amplitudes of 
both the maximum (at 5Tb ~ 38 mK) and minimum 
(at 5T b ~ 33 mK) of g2T>(5Tb) are higher than those of 
g2D(5Tb t n). Birth of new islands and peninsulas made of 
H I regions — some peninsulas may appear as islands at 
certain 5Tb — is responsible for the former, while birth 
of new lakes made of H II regions is responsible for the 
latter. Birth of islands and peninsulas occurs at mildly 
overdense regions (5Tb. n ~ 38 mK) because bubbles are 



clustered such that some of them merge with one another 
to fully or partly surround neutral regions. Note that a 
neutral region identified as an island in 2D may well be 
a cross-section of a vast neutral region in 3D. 

When the universe enters the middle phase of reion- 
ization (x v w 0.4 — 0.6), almost all overdense regions 
(5Tb, n ^ 30 mK) have been ionized to form larger bub- 
bles. This is reflected in the relatively low amplitude 
trough of g2v(5Tb), which appears in very low-density 
(5Tb j7l < 20 mK) regions. Now many more islands appear 
in mildly underdense (5Tb yU ~ 25 mK) regions because 
larger and more clustered bubbles penetrate further into 
low-density IGM and are more efficient in forming new 
islands. As time passes from x v = 0.4 to x v = 0.6, the 
amplitude of #2D decreases as bubbles merge with each 
other and neutral clumps disappear. Further penetra- 
tion of bubbles into lower-density IGM from x v = 0.4 to 
x v = 0.6 is also reflected in the maximum value of 5Tb 
for nonzero g2T>(5T\ ) ). 

Finally, in the late stage (x v = 0.9), all ionized bub- 
bles have been connected and the last surviving neutral 
clumps exist to give positive #2D at 5Tb ~ 20 mK. Note 
that these clumps exist only in underdense regions, which 
is a clear indication of the fact that reionization proceeds 
in an inside-out fashion: high-density regions will be ion- 
ized first due to the proximity to sources of radiation, and 
low-density regions will be ionized later. 

Now, we briefly investigate the impact of beam shapes. 
The impact of Wcg on #2D is demonstrated on the right 
panel of Figure [U Note that the Fourier transform of 
Wqg is proportional to k 2 cr 2 exp(— k 2 a 2 /2). From the 
convolution theorem, it is obvious that the filtered field 
becomes the Wc-smoothed, negative Laplacian (—V 2 ) 
of the field. Note that Laplacian is equivalent to the di- 
vergence of the gradient. Therefore, if the 5Tb field is 
filtered by Wqg, 5T b s in those regions with the steepest 
spatial gradient in x or 5 will correspond to extrema. We 
also note that the sign of the Laplacian depends on the 
morphology of x: if H II bubbles form in the sea of neu- 
tral gas, gradients of 5Tb diverge (— V 2 5Tb < 0), while if 
neutral clumps remain in the sea of ionized gas, gradients 
of 5Tb converge (— V 2 5Tb > 0). An interesting feature in 
g2T>(5Tf ) ) is indeed observed to move from regions with 
5Tb < at the relatively early (x v ~ 0.4) stage to regions 
with 5Tb > at the late (x v ~ 0.9) stage. Nevertheless, 
because the topology of the 21-cm signal is processed fur- 
ther (Gaussian smoothing and Laplacian) in the case of 
Wqg than that of Wg (Gaussian smoothing only), the in- 
terpretation becomes less transparent. We also note that 
Wcg filters out almost all scales except for the filtering 
scale, while Wq filters out only those scales smaller than 
the filtering scale. 

The nice interpretative power of g2T>(5Tb) is somewhat 
lost with Wcg- Based on this, we suggest that any ar- 
tificial effects from "dirty" beams should be removed or 
minimized in real observations. Only when the effec- 
tive beam shape becomes something close to Wq , general 
analyses including our genus method will obtain their full 
potential. 

4.3.3. Model- dependent feature and required sensitivity 

We finally investigate whether our 2D genus analysis 
can discriminate between different reionization scenar- 
ios. Ultimately, if this turns out true, we may be able to 
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Figure 6. Left: 2D genus of 5T b for fl25_1000S, filtered with 1' FWHM Gaussian beam and 0.2 MHz bandwidth (solid), compared to the 
template 2D genus out of the underlying density field (dashed), also filtered the same way. Right: Same as left, but with a compensated 
Gaussian beam applied. 
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probe properties of radiation sources at least indirectly, 
because these properties determine how cosmic reioniza- 
tion proceeds. 

Our 2D genus analysis seems very promising, indeed, 
in probing source properties. The top panel in figure [71 
shows gzD^Tb) and gzD{5Tb,n) of all cases at x v = 0.4. 
The impact of small- mass (10 8 < M/Mq < 10 9 ) halos is 
observed as follows. In fl25_1000S, almost all (filtered) 
overdense regions have been ionized at this epoch due 
to efficient, small-mass halos. In contrast, when there 
are no small- mass halos (f5, f250) or if small- mass halos 
are not as efficient (fl25_125S), some fraction of over- 
dense regions still remains neutral at x v ~ 0.4. These 
regions are mildly nonlinear, which are ionized due to 
the almost on-the-spot existence of efficient small-mass 
halos in fl25_1000S, while they mostly remain neutral or 
only partially ionized in f5, f250 and fl25_125S due to 
absence or low-efficiency of small-mass halos (figure fT0|) . 
The relatively larger amplitude of #2D in fl25_1000S is 
the reflection of the fact that there are many more small 
individual islands and bubbles, as explained in section 
14X21 

The impact of beam size and bandwidth is shown in 
figures [7] — As the applied Av increases, #2D curve 
shrinks in width and moves toward left. The interpre- 
tation we made on fl25_1000S with Av = 0.2 MHz and 
AO = 1' somewhat weakens in the sense that #2D curves 
are mapped too deep in (newly filtered) underdense re- 
gions. Nevertheless, the amplitude of #2D is the largest 
and the rightmost wing of #2D is located at the smallest 
5Tb in fl25_1000S in all varying filtering scales, just as 
when Av = 0.2'MHz is applied. 

Since #2D is found a useful tool for understanding the 
evolution (section [432]) and even discriminating between 
reionization scenarios, we need to ask whether this can be 
achieved in real observations. To estimate the required 
sensitivity, we show the degree of fluctuation in filtered 
5Tb in figure [TTJ 5Tb } rms becomes maximum around the 
middle stage of reionization in all cases. Roughly speak- 
ing, 5Tb, rms gives the required rms sensitivity limit of a 
radio antenna in each configuration, which is about a few 
mK. 

Following Ilie v et aTl (j2003), we estimate the sensitiv- 
ity limit of SKA for the filtering scale treated in this 
paper, to figure out feasibility of our analysis on the real 
data of 21-cm background. We assume the core size to 
be ~ 1 km for SKA. We also assume, just for the sake of 
sensitivity estimation, that x v ~ 0.5 is reached at z ~ 8. 
Even though our choice of input parameters (properties 
of sources) made reionization end much earlier than the 
usually believed epoch of end of reionization (z ~ 6.5), 
except for /5, we may imagine that similar distinctive 
features among different models would still exist if we 
tuned these parameters to make these models achieve 
the half- ionized state at z ~ 8. Figure fTTI shows the es- 
timated sensitivity under these conditions. Due to the 
strong dependence of sensitivity on ~ A# -2 , increase 
in the beam size quickly makes our analysis feasible. This 
comes at a price that values of g2F>(5Tb) shrinks further 
and the range of 5Tb shifts more to the left. At any rate, 
- 1000 - 10000 hours of integration at 5v ~ 1 - 2 MHz 
and 50 ~ 2' — 3' seems possible to allow our analysis. 

5. SUMMARY AND DISCUSSION 



We calculated the 21-cm radiation background from 
high redshift using a suite of structure formation and 
radiative transfer simulations with varying properties of 
sources of radiation. These properties of sources are spec- 
ified by the spectrum of halo masses capable of host- 
ing sources of radiation and the emissivity of hydrogen- 
ionizing photons. Assuming the pre-reionization heating 
of IGM, we calculated the 21-cm radiation background in 
the saturated limit, T$ ^> Xcmb, such that its differential 
brightness temperature is proportional to the underlying 
density and neutral fraction. 

In order to understand the topology of the high- 
redshift IGM, we developed a method of 2D genus topol- 
ogy applicable to the 21-cm radiation background we cal- 
culated. Basically, this method calculates the 2D genus, 
the difference between the number of hot spots and cold 
spots, under a given threshold differential brightness 
temperature. We construct 2D genus curves at different 
redshift s for different scenarios, by varying this threshold 
values. 

We found that the 2D genus curve gzD(5Tb), if com- 
pared to #2D of the underlying density fluctuation, 
#2D(^&,n)> clearly shows the evolution of the reioniza- 
tion process. g2r>(5Tb) is found to deviate quickly from 
the #2D of a Gaussian ra ndom field, qualita tively in ac- 
cord with the findings bv llliev et al . (2006) and m con- 
tradiction to the finding s by iShin. Trac fc Cenl (|2008D . 
iShin. Trac fc Cenl (|2008h find that the non-Gaussianity 
of the H II region is small even when the universe is 50% 
ionized. It is also shown that the reionization proceeds in 
an inside-out fashion: high-density regions gets ionized 
first and low-density regions gets ionized later. 

We also showed that our 2D genus method can be used 
to discriminate between various reionization scenarios, 
thus probing properties of sources indirectly. It seems 
most effective in discriminating the mass spectrums of 
halos which host sources of radiation. A hybrid mass 
spectrum with different emissivity in different species 
(fl25_1000S in this paper) stands out from cases with 
a single species (f5, f250) or a case with equal-emissivity 
among different species (fl25_125S). 

Crucial ingredients needed for this analysis is the beam 
shape and sensitivity. We tested two different beams, 
Gaussian and compensated Gaussian. The Gaussian 
beam, even after the field is filtered and degraded, leaves 
discernible imprint on the 2D genus curve such that dif- 
ferent reionization scenarios can be discriminated. In 
contrast, the compensated Gaussian beam somewhat 
loses the interpretative power that the Gaussian beam 
had. Assuming that the compensated Gaussian beam 
represents the so-called dirty beams, our analysis may 
be applied to full potential only when those artificial ef- 
fects from dirty beams are removed. 

We predict that SKA will be able to produce data suit- 
able for this analysis, when 1000 — 10000 hours of integra- 
tion is performed with 5v ~ 2 — 3 MHz and 50 ~ 2' — 3' 
at observing frequency of ~ 150 MHz. Therefore, our 
method seems promising. Note, however, that a direct 
link from the observed 2D genus curve to true proper- 
ties of sources will be still far-fetched, when we adopt 
such low-resolution filters to increase sensitivity. More- 
over, there still exist too many uncertainties, such as the 
matter power spectrum in small scales, which is rele- 
vant to formation and evolution of sources of reioniza- 
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Figure 7. 2D genus curves of the differential brightness temperatures at x v = 0.4, filtered with V FWHM Gaussian beam and with 
0.2 MHz (top), 1 MHz (middle), 2 MHz (bottom) bandwidth, compared to those of the matter densities (black). 



tion. We hope to see many more useful constraints on 
cosmic reionization come from various other, direct or in- 
direct, observations. In the future, we will explore more 
reionization scenarios to strengthen the potential of our 
2D genus analysis. 
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Figure 9. Same as figure [7] but with 3 ; FWHM Gaussian beam applied. 
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Figure 10. How 2D genus depends on reionization scenarios: from left to right, threshold ST^s (contour lines) are selected which give the 
minimum ^D^T^), maximum <72d(<^&) 5 an d zero <72D(<^6,n)j respectively, together with the snapshot of unfiltered ionized fraction (far 
right). x v = 0.4 in both (top: fl25_1000S; bottom: fl25_125S) cases. See arrows in figure [3 for contours a, b and c. 




Figure 11. Top row: Sensitivity limit of SKA in terms of bandwidth Au, with varying integration time (100, 1000 and 10000 hours). Also 
plotted are the required sensitivity limits (cross) for useful results, at Au = 0.2, 1 and 2 MHz. Each column corresponds to varying beam 
size A#'s (top labels). The rest: Evolution of rms differential brightness temperature for varying AO (each column) and Au (each row). 



